Can half metallic zincblende MnAs be grown? 
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In this paper we investigate theoretically the ground state NiAs-type structure of MnAs and 
we compare the magnetic and structural properties with an hypothetical zincblende structure. 
A zincblende structure can be obtained, in principle, from the diluted magnetic semiconductor 
Gai-aiMua; As in the limit x = 1. Using density functional calculations within the local spin-density 
approximation (LSD A), we show that the zincblende structure, although showing half metallic be- 
havior which is very attractive for spintronics, can not be stabilized at equilibrium. We perform a 
tight-binding analysis of the Mn-As bond in the tetrahedral coordination to investigate the nature 
of the bonding in diluted magnetic semiconductors. 
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PACS: 71.15.H, 71.20.M, 75.50.P, 71.15.F 



I. INTRODUCTION 



In the last few years there has been a growing interest 
in combining magnetotransport experiments based onto 
the giant magnetoresistance (GMR) in metallic magnetic 
multilayers with semiconductor physics Semicon- 
ductors offer various advantages over metals such as a 
very long spin- hfetime |], persistent spin coherence 
and compatibility with existing semiconductor process- 
ing technologies. These properties pave the way for us- 
ing semiconductors as media for storing coherence and 
realizing elements for solid-state quantum computation 
Q . The natural extension of GMR to semiconductors is 
in a hybrid spin-valve where the magnetic elements are 
either magnetic metals or diluted magnetic semiconduc- 
tors IJI and the non-magnetic elements are ordinary 
semiconductors. In both cases it is essential to inject 
spins into the non-magnetic semiconductors and this is a 
formidable challenge. 

So far spin-injection into semiconductors from metallic 
contacts has been elusive [Q with GMR signals smaller 
than 1% In contrast spin-injection from diluted 

magnetic semiconductors has been more successful and 
evidence of polarized currents either into III-V and 
II- VI semiconductors has been demonstrated. In 
both cases the polarization is observed by optical tech- 
niques, which makes the integration of such systems in 
present semiconductor technology very difficult and all 
electronic systems would be more technologically desir- 
able. Recently it has been suggested jlj] that a spin- valve 
with a conventional injector-detector geometry made by 
metallic contacts cannot present large signals at equilib- 
rium due to the very different resistivities of the contacts 
and the semiconductor spacer. Much larger signals are 



predicted by using half-metallic contacts or diluted mag- 
netic semiconductors. 

From the point of view of the materials the system 
formed by GaAs as non-magnetic semiconductor and by 
MnAs and Gai_a;Mna;As respectively as magnetic metal 
and diluted magnetic semiconductor is the most promis- 
ing for spintronics applications. For the metallic compo- 
nent in metal/semiconductor structures, MnAs presents 
several advantages over transition metals. MnAs grows 
epitaxially onto GaAs in the hexagonal NiAs-type struc- 
ture. This structure is characteristic of all bulk MnX 
compounds with X=As, Sb and Bi. High quality all- 
epitaxial MnAs structures can be grown onto commonly 
[l5| and Si Q. For MnAs/GaAs the in- 



used GaAs [13 



terfaces are thermodynamically stable since the materials 
share the As atoms and the growth process is completely 
compatible with existing III-V MBE technology. 

Another attractive prospective crystal structure for 
MnAs is the zincblende structure. Recent density func- 
tional calculations have shown that the latter, which 
can be obtained in principle in the limit of 100% dop- 
ing from Gaa;Mni_a;As, is a half-metal. This makes 
zincblende MnAs very promising for spintronics because 
it will allow the intrinsic difhculty of injecting spins into 
semiconductors to be overcome [|l2|. Unfortunately to 
our knowledge zincblende MnAs has never been success- 
fully grown. 

In this paper we investigate, by using density func- 
tional calculations, the stability of the NiAs-type MnAs 
with respect to the hypothetical zincblende structure. In 
particular we look for growth conditions (for example 
strains from the growing substrate, volume compressions) 
which can enable the zincblende phase to be made. One 
of the main results of this analysis is that the zincblende 
structure cannot be stabilized either by lattice stretch- 
ing nor by compression. We also consider the magnetic 
properties of the NiAs-type MnAs grown onto GaAs with 
different crystalline orientations and lattice distortions. 
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Finally we look at the nature of the chemical bonding in 
zincblende MnAs. This is particularly important since 
the tetrahedral coordination of the Mn atoms is one of 
the key elements for understanding the properties of mag- 
netic diluted III-V semiconductors. 

The remainder of this paper is organized as follows. In 
the next sections we will briefly discuss the calculation 
technique and provide some information regarding the 
lattice structures considered and their properties. Then 
we will move on to analyze the stability of the NiAs- 
type MnAs under compression and distortion of the unit 
cell. In section V we will consider the zincblende struc- 
ture and we will compare its structural properties with 
the NiAs-type. Finally we will perform a fit onto a tight- 
binding model in order to analyze the details of the chem- 
ical bonding in zincblende MnAs. 



II. COMPUTATIONAL DETAILS 

We compute the electronic structure of MnAs with 
both NiAs-type and zincblende lattice structures using 
a plane wave pseudopotcntial implementation of density 
functional theory ||l8| with the local spin density approx- 
imation. The use of pseudopotentials for the study of 
magnetic systems, although not yet as popular as all 
electron methods, is now well documented We use 
the optimized pseudopotentials developed by Rappe et 
al. for our Mn pseudopotcntial, allowing us to use a 
small cutoff energy for the plane wave expansion in the 
total energy calculation of the solid. This is achieved by 
minimizing the kinetic energy in the high Fourier com- 
ponents of the pseudo wave function. The reference con- 
figuration for our Mn is 4s*'-^^4p'''^^3d^ with core radii of 
2.0, 2.15 and 2.0 a.u. and we use the usual Kleinman- 
Bylander separable form with two projectors for each an- 
gular momentum |^,^. For As we use the standard 
Hamann-Schiiltcr-Chiang pseudopotcntial [p3|| with one 
projector for each angular momentum. The pseudopo- 
tentials were tested by comparing the pseudoeigenvalues 
with those generated by all electron calculations for sev- 
eral atomic and ionic configurations. We also checked 
the portability by calculating the equilibrium lattice con- 
stant of GaAs and the dependence of the magnetization 
on the unit cell volume of hypothetical fee Mn. Both 
are in excellent agreement with all electron calculations. 
Moreover the Mn pseudopotentials has been successfully 
used in earlier studies of perovskite manganites | p9[ . 

Total energy and band structure are calculated by us- 
ing the code SPECTER [||. This is a spin-polarized 
implementation of density functional theory written in 
C, which originated from the program CASTEP 2.1 [^ . 
The plane wave cutoff is fixed to 870 eV (~63 Ry) and we 
use a simple electron density mixing scheme for conver- 
gence developed by Kerker with variable Gaussian 
broadening between 0.8 and 10~^ eV. We use a 4 x 4 x 4 



Monkhorst-Pack grid for the zincblende structure and a 
3 X 3 X 3 for the NiAs-type structure. This leads re- 
spectively to 10 and 6 fc-points in the corresponding ir- 
reducible Brillouin zones. Stability of the results with 
respect to fc-point sampling has been carefully checked 
and the present choice represents a good compromise be- 
tween accuracy and computation time. 



III. CRYSTAL STRUCTURES 

The NiAs-type lattice structure is an hexagonal struc- 
ture (space group PBs/mmc) with 4 atoms in the prim- 
itive unit cell (see figure It is the equilibrium lattice 
structure of the pnictides formed with Mn at room tem- 
perature. Bulk MnAs has lattice constants of a=3.7A 
and c=5.7A and is a ferromagnetic metal with a Tc of 
318K. Above Tc it undergoes a phase transition to a 
paramagnetic state involving also a structural change to 
the hexagonal MnP-type lattice structure (space group 
Pnma) with a 2% volume reduction |27|. A further struc- 
tural phase transition back to the NiAs-type phase is 
observed at 398K with the material remaining paramag- 
netic. 
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FIG. 1. NiAs-type crystal structure. The lattice constants 
for bulk MnAs are a=3.7A and c=5.7A . 



MnAs can grow epitaxially onto GaAs in sev- 
eral possible crystalline orientations with respect to 
the GaAs substrate, depending on the growth con- 
ditions. In particular when MnAs grows onto (001) 
GaAs two orientations can occur The first one, 

which is usually called type A, has growth plane 
(1100) and epitaxial relationship [1120]MnAs|| [110]GaAs 
and [0001]MnAs||[110]GaAs, while the second, known 
as type B, has growth planes (1101) and (1102) 
and epitaxial relationship [ii20]MnAs|| [llOjGaAs and 
[li20]MnAs|| [llOjGaAs. Moreover MnAs can be grown 
onto (111) GaAs ||2^] with growth direction along the 
c axis of the hexagonal cell and epitaxial relation- 
ship (0001)MnAs||(lll)GaAs, [1100]MnAs|| [112]GaAs 
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and [2110]MnAs|l[110]GaAs. All these phases have the 
easy axis of magnetization along the c axis of the hexag- 
onal cell. The saturation magnetizations (M^) are very 
different, being 416 emu/cm^, 331 emu/cm'^ and 
640 emu/cm'^ respectively for type A, type B and 
(111) GaAs grown MnAs. Moreover for MnAs MBE- 
deposited onto GaAs the saturation magnetization de- 
pends on the number of monolayers deposited ||30|] , sug- 
gesting that relaxation of the lattice has an important 
influence on the magnetic properties. In what follows we 
will show that these changes of the saturation magnetiza- 
tion can be correlated with distortions of the hexagonal 
unit cell resulting from matching to the substrate. 

Bulk zincblende MnAs is unstable, and zincblende 
structure thin films have not been demonstrated con- 
clusively |^l|. However Mn can be incorporated into 
epitaxial layers of GaAs with concentrations beyond its 
solubility limit by using low temperature MBE tech- 
niques The highest concentration obtained so far 
is cc = 0.07 for Gai_2:Mn2:As, above which the Mn atoms 
segregate forming clusters of MnAs with Ni As- type struc- 
ture. Nevertheless the low concentration limit is repro- 
ducible. Mn in GaAs acts both as a source of localized 
spin and also as an acceptor, providing electrical carriers 
(holes). (Ga,Mn)As is ferromagnetic for concentration 
above x = 0.005 although the origin of the ferromag- 
netism is still a matter of debate. In section V we will 
study the limit of x = 1 in order to understand the nature 
of the chemical bond of Mn tetrahedrally coordinated 
with As. 



IV. NIAS-TYPE MNAS 



In this section we look at the magnetic properties and 
structural stability of the NiAs-type phase. In figure || 
we present the band structure for ferromagnetic MnAs 
with lattice parameter corresponding to the experimen- 
tal bulk values (a=3.7A c=5.7A ). The two lower bands 
correspond to the 4s states of As and the following 16 
to hybridized p-d bands. This band structure is char- 
acteristic of all the Mn-based pnictides and we can 
identify three distinct regions: i) a low energy p-d bond- 
ing region (between -6 and -2.5 eV for the majority band 
and between -6 and eV for the minority), ii) an almost 
dispersionless intermediate region (~-2.0 for the majority 
and ~ 1 for the majority band) and iii) a high energy p-d 
antibonding region (between -1.5 and 3 eV for the ma- 
jority band and between 1.5 and 4 eV for the minority). 
The ferromagnetism is mostly due to the spin-splitting 
of the degenerate d bands, although a rigid shift model 
cannot be applied due to the strong p-d interaction. We 
will see in the next section that p-d interaction is also 
very strong in the zincblende structure. 
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FIG. 2. Band structure for bulk MnAs with NiAs-type lat- 
tice structure. The figure on the left corresponds to the ma- 
jority spin and the one to the right to the minority. The hor- 
izontal line denotes the position of the Fermi energy, which 
has been chosen to be eV. 

Turning our attention to the structural properties we 
calculate the total energy and the magnetic moment per 
MnAs pair as a function of the volume occupied by a pair 
of atoms. We fix the ratio between the hexagonal axes to 
c/a=1.54, which corresponds to the experimental value 
for bulk MnAs. The results are presented in figure 0. 
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FIG. 3. Total energy (□ on the left-hand side scale) and 
magnetization ( x on the right-hand side scale) as a function 
of the volume occupied by a MnAs pair. 

Our calculated equilibrium volume is 31.123 A'^ which 
is slightly smaller than the experimental one of 33.789 A3. 
This small discrepancy can be explained by a general 
tendency of the LDA to underestimate the equilibrium 
volume. We also note that the structure is very sta- 
ble with respect to compression and expansion of the 
unit cell around the equilibrium value (volume changes 
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of the order of 50% give energy changes of approxi- 
mately 1 eV). Our calculated magnetization Mg at equi- 
librium is between 2.5 and 3.0/LtB which is smaller than 
that found experimentally (~ 3.4/LtB). Nevertheless it is 
notable that the magnetization is very sensitive to the 
volume around the equilibrium position, and small vol- 
ume changes can produce large changes in the magneti- 
zation. For large compressions the system evolves to a 
paramagnetic state. This transition is not directly com- 
parable with experiments on MnAs under pressure js^ j 
since a structural transition to the MnP-type structure 
and a magnetic transition to a mixed ferromagnetic and 
antiferromagnetic phase are experimentally observed for 
pressures above 4 kbar. In contrast, for volumes larger 
than the equilibrium volume we observe a saturation of 
the magnetic moment per MnAs pair to a value around 
4//B. This saturation can be understood by comparing 
the band structures of Fig.0 with those obtained for an 
expanded structure (a=4.3A c=6.62A with a volume per 
MnAs pair of 53.03A'^) presented in Figj^. First we ob- 
serve that there is a large shrinking of the Mn d bands 
without a drastic change of their centers with respect to 
the Fermi energy. As a result of the strong p-d inter- 
action the minority p bands below the Fermi energy are 
pushed to higher energies, while the majority p bands 
above the Fermi energy are pulled to lower energy. This 
increases the number of occupied majority states and de- 
creases the occupation of minority states resulting in a 
global increase in the magnetic moment of the system. 
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FIG. 4. Band structure of NiAs-type MnAs at expanded 
volume. The lattice constants are a=4. 3A and c=6.62A . The 
figure on the left shows the majority spin and the one to the 
right the minority. The horizontal line denotes the position 
of the Fermi energy, which has been chosen to be eV. 



The saturation of the magnetization to a value of 
~ 4/iB can be explained by considering charge transfer 
from Mn to As. Three electrons will be transferred from 



each Mn to a neighboring As in order to fill completely 
the electronegative As atoms' 4p shell. This leads to an 
atomic configuration of Mn with four electrons in the va- 
lence, which arrange themself following Hund's rules and 
produce a magnetic moment of A^b- The charge trans- 
fer is not complete for unit cell volume corresponding to 
the experimental value because of the strong p-d inter- 
action. This p-d interaction in fact leads also to a quite 
large negative polarization of the anions (-0.23±0.05^b), 
which has been observed with neutron scattering [Q. 
For stretched unit cells, the p-d hybridation is reduced 
and a more ionic configuration is resumed. 

We finally look at the magnetization and the total en- 
ergy of distorted cells. We consider cells with a volume 
equal to that of bulk MnAs and calculate the energy and 
the magnetization as a function of the cell aspect ratio 
c/a. Our results are presented in Fig.^j 
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FIG. 5. Total energy (□ on the left-hand side scale) and 
magnetization ( x on the right-hand side scale) as a function 
of the cell aspect ratio c/a. For bulk MnAs c/a=1.54. 

Several interesting aspects must be pointed out. First 
note that the energy scale is much smaller than that of 
Fig.^. The total energy for aspect ratios from 1.4 to 1.7 
changes only by 0.25 meV (which is equal to fc^T at 
room temperature). By comparing figures I and I one 
can conclude that NiAs-type MnAs can accommodate 
large distortions, which are much less energy demanding 
than volume changes. This can be easily understood by 
using a simple geometrical argument. Consider the near- 
est neighbor distance dMn-As as a function of the lattice 
constant a for fixed cell volume 



-^Mn-As 



(1) 



VsV^ a6' 

where Vb is the volume occupied by a Mn-As pair. It 
clear that dMn-As is non-monotonic with a minimum at 
a = (2Vb)^/^, which is 3.629A in the case of bulk MnAs. 
Moreover dMn-As docs not vary much with the lattice 
constant a around this minimum. For c/a increasing 
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from 1.3 to 1.8, a increases by about 10% while dMn-As 
increases by < 1%. Since the total energy does not vary 
considerably as a function of the Mn-As bond angle but 
only as a function of its bond length, we can conclude 
that the Mn-As bond is largely ionic and dominates the 
total energy. In fact for the lattice distortions considered 
the second nearest neighbor distances (As- As and Mn- 
Mn) change drastically, while the first nearest neighbor 
distance (Mn-As) is almost constant. This produces no 
appreciable change in the total energy suggesting that the 
latter is dominated only by the Mn-As bond. Moreover 
because of the absence of any strong angular dependence 
we can also conclude that the Mn-As has a large ionic 
component. 

Our calculated tiny energy variation for even large dis- 
tortions, combined with the large change in magnetiza- 
tion at small distortions explain most of the properties of 
the different epitaxial orientations of MnAs onto GaAs. 
In fact it is worth noting that a large change of the mag- 
netization occurs for ratios c/a just above the one cor- 
responding to the bulk. The magnetization passes from 
a high magnetization region c/a<1.54 to a low magne- 
tization region c/a>1.54. This is very important since 
large changes of the magnetization have been observed in 
MnAs grown onto GaAs with different crystalline orien- 
tations (type A, type B and MnAs onto (111) GaAs). It 
is difficult to correlate directly the deformation of the cell 
with the magnetization observed in experiments since a 
complicated surface reconstruction occurs and a de- 
tailed characterization is not easy. We do not want to go 
into details of this intricate issue and we only suggest that 
the reduction of the magnetization of MnAs grown onto 
GaAs may be due to large deformation of the hexagonal 
lattice. 



V. ZINCBLENDE MNAS 

We now consider MnAs with the zincblende lattice 
structure. In Fig.||we present the total energy as a func- 
tion of the lattice constant for MnAs and also compare 
the total energy of the same structure when the mag- 
netic moment is fixed either to O/i^ (paramagnetic phase) 
and to 4/i b (half- metallic phase) . The equilibrium lattice 
constant is found between 5.6A and 5.7A which is very 
similar to the lattice constant of GaAs. Such a value is 
smaller (as we would expect within the LDA) than that 
predicted from the linear extrapolation of the experimen- 
tal lattice constant of Gai_a;Mn2:As for cc — > 1 (5.89A) 
[^. It is also smaller than another LDA result obtained 
with the FLAPW method (5.9A) H]. This could be due 
to the pseudopotential method used here or differences 
in convergence or sampling between the two calculations. 
We also notice that the paramagnetic phase becomes en- 
ergetically favorable for a=5.0A, where a ferromagnetic 
to paramagnetic transition is predicted. In contrast, for 



lattice spacings larger than 5.8A a transition to an half- 
metallic state is predicted. This is a very important result 
since half-metallicity is one of the conditions required to 
obtain large signals in a semiconductor spin- valve [E2I. 
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FIG. 6. Total energy as a function of the zincblende lat- 
tice spacing for the paramagnetic phase (•) , the half- metallic 
phase (*) and for non-constrained magnetic moment (□). 

The behavior of the magnetization as a function of the 
unit cell volume is summarized in Fig.|^. 
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FIG. 7. Total energy (□ on the left-hand side scale) and 
magnetization ( x on the right-hand side scale) as a function 
the unit cell volume. The half-metallic behavior occurs for 
a > 5.8A. 

We notice that Fig.|^ looks similar to its counterpart for 
the NiAs-type structure (Fig.^) although in the latter no 
transition to an half-metallic state is found. This funda- 
mental difference between the two structures can be un- 
derstood by looking at the band structure of zincblende 
MnAs in Fig.|. 

We first note that, excluding the presence of the Mn 
d bands, the bandstructure closely resembles that of the 
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non-magnetic III-V semiconductors. Consider the ma- 
jority bands first. By symmetry analysis we can easily 
identify the lower lying As s states (lowest band, Fi at 
the r point), the As p valence band (first Fis point above 
Ef) and the first of the conduction bands (first Fi point 
above Ep). However, due to the strong interaction with 
the Mn d states, the As p bands (top of the valence bands 
in GaAs) are pushed toward higher energies and become 
half filled. The Mn d bands, which are split into the dou- 
bly degenerate Cg band (F12) and the triply degenerate 
t2g band (F15) are below the Fermi energy and entirely 
occupied. Therefore one can conclude that two impor- 
tant interactions are present: i) the s-p interaction giving 
rise to the sp^ bond typical of non-magnetic semiconduc- 
tors and ii) the p-d interaction responsible for the mag- 
netism, which sets the position and the dispersion of the 
As p band. Note also that in contrast to the prediction 
from crystal field theory, the eg orbitals possess higher 
energy than the t2g. This is also due to the strong p-d 
interaction. If we now turn our attention to the minority 
band, we can find most of the features of the majority. 
The main difference is the large split between the t2g and 
the Cg states, which gives rise to a large gap in the band 
structure. 




-12 I 

Lr XK TLT XK T 

FIG. 8. Band structure for zincblende MnAs at the LDA 
energy minimum (a=5.7 A). Thefi gure on the left corresponds 
to the majority spin and the one on the right to the minority. 
The horizontal line denotes the position of the Fermi energy, 
which has been chosen to be eV. 



We observe that the Fermi energy cuts through the 
band center of the high dispersion p band in the ma- 
jority band, and through the band edge of the almost 
dispersionless d band in the minority. This is very impor- 
tant from the point of view of the transport, since very 
different spin-dependent effective masses are predicted. 
When the lattice constant is increased, we expect a gen- 



eral narrowing of all the bands. This does not affect 
qualitatively the majority band. However for some crit- 
ical lattice spacing the band edge of the Cg band in the 
minority spin shifts above Ep giving rise to a semicon- 
ducting behavior. This generates the half metallicity. 

In summary, if we start from the atomic configura- 
tion for Mn and As, the mechanism giving rise to the 
magnetism is the transfer of one electron from the Mn d 
orbitals to the As p shell in order to form the sp'^ bond. 
The remaining 4 electrons in the Mn d shell maximize 
the local magnetic moment of Mn due to Hund's cou- 
pling and the final state turns out to be ferromagnetic 
due to the strong p-d interaction. Note that these are 
features of the tetrahedral coordination and are absent 
in the NiAs-type structure, where no half metallic state 
is predicted. 

Finally we look at the stability of the zincblende struc- 
ture with respect to the NiAs-type. The total energy and 
the saturation magnetization of both the structures as a 
function of the volume occupied by a MnAs pair are pre- 
sented in Figjsl. 
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FIG. 9. Total energy (a) and magnetization (b) as a func- 
tion of the MnAs pair volume for the zincblende (•) and 
the NiAs-type structure (□). Note the large stability of the 
NiAs-type structure over a very broad volume range. 

It is clear that the NiAs-type has a much lower total 
energy and also a more compressed lattice. Therefore it 
is the stable structure at all the thermodynamically ac- 
cessible pressures. This is consistent with the very small 
dilution limit of Mn in zincblende GaAs and with the fact 
that the annealing usually induces segregation of MnAs 
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NiAs-type particles within the GaAs matrix ||^,^. 

We also note that for very expanded unit cells the 
zincblende structure becomes stable with a crossover vol- 
ume which corresponds to a zincblende lattice constant 
of 5.8A. However this large volume increase in the NiAs- 
type structure is very unlikely to occur, even if the lattice 
constants are forced to be large by growing onto sub- 
strates with large mismatch. As we noted previously, the 
NiAs-type structure can easily accommodate large cell 
distortions (Fig.^. Therefore, at equilibrium it is ener- 
getically more favorable for the system to distort the cell, 
instead of increasing the volume and inducing a NiAs- 
type to zincblende transition. 

Finally we investigate the effect of the anion size on 
the stability properties of the zincblende structure. Our 
hope is that for larger ionic radii the NiAs-type struc- 
ture should become more rigid and a transition to the 
zincblende structure at large volume may be induced. 
We perform the same calculations described above for 
MnBi both with NiAs-type and zincblende structure 
(bulk MnBi is a ferromagnetic metal at room temper- 
ature with NiAs-type structure and lattice constants 
a=4.17A, c=5.764A). We did not find any relevant qual- 
itative difference with respect to the MnAs case (see 
Fig.|l^) and this is related to the geometrical argument 
presented in section IV. Hence we conclude that, for all 
the Mn-pnictides, the NiAs-type is the more stable crys- 
talline structure and a transition to zincblende structure 
can be obtained only with highly non-equilibrium meth- 
ods. 
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FIG. 10. Total energy as a function of the MnBi pair vol- 
ume for the zincblende (•) and the NiAs-type structure (□). 
Note the large stability of the NiAs-type structure over a very 
broad volume range as in the MnAs case. In the inset the 
energies for three different values of c/a for a volume cor- 
responding to the bulk MnBi are shown. The central point 
corresponds to bulk MnBi {V = 43.37A^). 



VI. MN-AS BOND IN THE ZINCBLENDE 
STRUCTURE 



Finally, to analyze the detailed nature of the Mn-As 
bonding in the zincblende lattice we perform a tight- 
binding fit to our calculated LDA band structure. We 
use a fitting algorithm included in the package OXON 
p6| which minimizes the following function 



(2) 



where 7 is the m-dimensional vector containing the 
Slater-Koster tight-binding parameters a„ is the 

weight assigned to the eigenvalue En and E^ is the com- 
puted eigenvalue. We consider first (Mn-As) and second 
(Mn-Mn and As- As) nearest neighbor couplings. This in- 
teraction set has been shown to reproduce correctly the 
valence band and the first of the conduction bands of 
diamond-like semiconductors |^^. In the case of MnAs 
the presence of the Mn d band makes the fitting pro- 
cedure more complicated then in non-magnetic semicon- 
ductors. However we do not need to introduce any ad 
hoc excited state (such as s*) to fit the valence band 
and the first conduction band as is frequently done in 
the literature . The fit has been performed using 350 
eigenvalues which correspond to the lowest 10 energies 
calculated at 35 fc-vectors and the weight is chosen to 
be 1 for each eigenvalue. We perform two independent 
fits for the majority and minority band allowing both the 
on-site energies and the hopping integrals to be different. 

In figures |Tl| and |l2| we present the best fits for the 
majority and minority spin band respectively. We note 
that the agreement with the ab initio bands is remarkably 
good, in particular for the lower lying bands. The high- 
est bands are not well reproduced (in particular for the 
majority spin) because we did not include higher energy 
bands in the fit. The complete parameterization giving 
rise to the bands in figures |l^ and ^ is provided in tables 
|and|l| 

Several important points must be stressed. First we 
note that the major difference between the majority and 
the minority parameters is the value of the Mn 3d and As 
4p on-site energies. The on-site energy of the 3d Mn or- 
bitals coincides with the position of the eg states, which 
by symmetry are weakly coupled to the other bands. By 
contrast the t2g states are strongly coupled to the above 
As p band and this coupling is strongly spin dependent. 
From table || one can see that the dominant hopping 
integrals are along the Mn-As bond and that the param- 
eters are quite similar for the two spin directions. The 
only parameters which differ strongly for up and down 
spin are the sdcr, pdcr and pdvr elements. These govern 
the repulsion of the t2g Mn d bands and the As p bands 



7 



and are responsible for the large eg-t2g splitting in the 
minority band. 




L r X u,K r 

FIG. 11. Fit of the band structure for the majority band 
of MnAs. The stars (*) are the ab initio eigenvalues and the 
continuous line are the band calculated with the parameters 
of Tables | and [| 




L r X u,K r 

FIG. 12. Fit of the band structure for the minority band 
of MnAs. The stars (*) are the ab initio eigenvalues and the 
continuous line are the band calculated with the parameters 
of Tables § and [| 

Finally we check our hypothesis on the formation of a 
strong sp^ bond by comparing the ss, sp and pp hopping 
integrals for MnAs with those of GaAs, which have been 
obtained with the same fitting procedure from our cal- 
culated LDA band structure. The GaAs parameters are 
shown in table III and their magnitude is clearly com- 



parable with the MnAs case. Moreover also the second 
nearest neighbor pptr and pp7r integrals, which are the 
only second nearest neighbor parameters with apprecia- 
ble amplitude and the ones that describe most of the 
dispersion of the As p band along K— > F, are very sim- 
ilar. This clearly suggests the formation of a strong sp^ 
bond in MnAs together with a strong p-d hybridation. 



VII. CONCLUSION 

We have investigated theoretically the structural and 
magnetic properties of the hexagonal NiAs-type and the 
zincblende phases of MnAs. From the analysis we con- 
clude that the NiAs-type structure is more stable at all 
thermodynamic pressures. A phase transition to the 
zincblende structure is predicted for volume stretching, 
which in principle could be obtained by growing onto a 
substrate with large lattice mismatch. However the NiAs 
phase can accommodate large lattice distortions which 
preserve the Mn-As nearest neighbor distance and a tran- 
sition to a zincblende phase looks very unlikely. 

We have also investigated the zincblende phase in more 
detail, since the tetrahedral coordination is crucial for the 
physics of diluted magnetic semiconductors. By perform- 
ing both ab initio calculations and tight-binding fitting 
we see that two important effects are present: i) the for- 
mation of a stable sp^ bond and ii) a strong p-d hybri- 
dation, which splits the Mn d bands and gives rise to 
ferromagnetism. In particular for large lattice spacings 
the zincblende structure is predicted to be half metallic. 
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Majority (eV) Minority (eV) 



Eas{4s) -9.8394 -9.9188 
Eas{4p) -0.8709 -0.3476 
EMn{4s) 0.7380 -1.0519 
EMn(4p) 1.1254 1.4601 
EMn(3d) -5.4699 -2.3417 



TABLE I. On site energies for the majority and minority 
bands of zincblende MnAs. The notation is that of Slater and 
Koster [37]. 







Majority (eV) 


Minority (eV) 


E[As(4s) 


— Mn(4s)]£j 


1.8554 


1.3592 


E[As(4s) 


-Mn(4p)]o- 


-1.8684 


-1.5061 


E[As{4s) 


-Mn(3d)]CT 


1.0852 


2.0521 


E[As(4p) 


— Mn(4s)]CT- 


2.5898 


2.3806 


E[As(4p) 


— Mn(4p)]CT 


2.2825 


2.5146 


E[As(4p) 


— Mn(4p)]7r 


-0.9342 


-0.8346 


E[A3{4p) 


-Mn(3d)]<T 


1.0329 


1.4592 


E[As(4p) 


-Mn(3d)]7r 


-0.4317 


-0.6831 


E[As(4s) 


— As(4s)]cr 


-0.0133 


-0.1377 


E[As{4s) 


— As(4p)]o- 


0.0031 


0.0175 


E[As(4p 


-As{4p)]o- 


0.0009 


0.0064 


E[As(4p 


— As(4p)]7r 


0.0255 


-0.0373 


E[Mn(4s) 


— Mn(4s)]o- 


0.0000 


0.0000 


E[Mn(4s) 


— Mn(4p)]CT 


0.0053 


-0.0009 


E[Mn(4s) 


-Mn(3d)]<T 


-0.1579 


0.1216 


E[Mn(4p) 


— Mn(4p)]CT 


0.4411 


0.3323 


E[Mn(4p) 


— Mn(4p)]7r 


-0.2265 


-0.0167 


E[Mn(4p) 


-Mn(3d)]a 


0.2398 


-0.2232 


E[Mn(4p 


-Mn(3d)lir 


-0.0298 


0.1963 


E[Mn(3d 


-Mn(3d)]<T 


-0.1139 


-0.1199 


E[Mn(3d 


-Mn(3d)]7r 


0.0453 


0.0773 


E[Mn(3d 


-Mn(3d)l(5 


-0.0002 


-0.0005 





TABLE II. Hopping integrals for the majority and minor- 



ity bands of zincblende MnAs. The notation is that of Slater 
and Koster [37]. 
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GaAs (eV) 



Eas{4s) 


-7.7859 


Eas(4p) 


1.0395 


Ec;a(4s) 


-0.2676 


EGa(4p) 


2.8029 


E[Aa(4s)-Ga(4s)]CT 


1.4895 


E[As(4s)-Ga(4p)]a 


-1.5445 


i!^[As(4p)-Ga(4s)]cr 




E[As(4p) — Ga(4p)](7 


2.1571 


E[As(4p)-Ga(4p)]7r 


-0.7540 


E[As{4s)-As(4s)]ct 


0.0032 


E[As(4b)-As(4p)]ct 


-0.0437 


E[As(4p)-As(4p)]CT 


-0.0009 


E[As{4p)-As(4p)]7r 


0.0010 


E[Ga(4s)-Ga(4s)l<T 


-0.1038 


E[Ga(4s)-Ga(4p)]CT 


-0.0919 


E[Ga(4p)-Ga(4p)]o- 


0.4828 


E[Ga(4p)-Ga(4p)]7r 


-0.1512 


TABLE III. Tight-binding parameters 


for GaAs. The no- 


tation is that of Slater and Koster [37] . 





